Self-bound models of compact stars and recent mass-radius measurements 
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The exact composition of a specific class of compact stars, historically referred to as "neutron 
stars", is still quite unknown. Possibilities ranging from hadronic to quark degrees of freedom, 
including self-bound versions of the latter have been proposed. We specifically address the suitability 
of strange star models (including pairing interactions) in this work, in the light of new measurements 
available for four compact stars. The analysis shows that these data might be explained by such an 
exotic equation of state, actually selecting a small window in parameter space, but still new precise 
measurements and also further theoretical developments are needed to settle the subject. 

PACS numbers: 97.60.Jd, 26.60.Kp, 21.65. Qr 
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I. INTRODUCTION 

The first reasonable ideas about tiie composition of 
compact stars stated that matter under extreme densi- 
ties is composed of hadrons, i.e., mainly neutrons with 
small fractions of protons and also electrons. Further the- 
oretical developments and modern experimental results 
opened the window to other possibilities as the existence 
of hyperons inside such stars due to the high densities, 
and a whole new class of condensates (see, for example, 
[l| and references therein). Another possibility arisen in 
the '70s is the existence of deconfined quarks inside com- 
pact objects, either located only in the inner core of these 
stars or present up to their surfaces. The latter extreme 
possibility was suggested about three decades ago 0-[ai) 
emerging as an astrophysical realization of the stability 
scenario of the so-called strange matter. 

Consequences of the so-called strange matter hypoth- 
esis for compact stars have been extensively analyzed 
within the simple framework of the MIT Bag Model [3, Q 
and the Nambu- Jona-Lasinio model |9l-[ll| (see [l^l for an 
overview on the subject). After an initial round of per- 
turbative considerations [Tj, [lj| , the possibility of non- 
perturbative pairing between quarks [^, Il5l - [l9| allowed a 
new view on this subject due to a great enhancement on 
the window of stability for strange quark matter (having 
the strange quark mass, bag constant and the novel pair- 
ing energy of the quark condensate as parameters in this 
approach) . This state of totally paired 3-flavor quarks is 
called the color-flavor locked (CFL) strange quark mat- 
ter, in which quarks form Cooper pairs, thus lowering 
further the energy of the system. 
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Much work has been put forward in order to have a bet- 
ter understanding on the actual composition of neutron 
stars, both observational and theoretical. Nevertheless, 
there is no absolute conclusion for any model as yet. 

We intend in this work to contribute to the understand- 
ing of how the new and much more precise astrophysical 
measurements of the mass and radius of neutron stars 
[2(]| - l23| can help revealing the viability of exotic quark 
star models. We analyze in particular the viability of 
the class of CFL strange quark matter models. They are 
likely to be the most favorable candidates for self-bound 
stars, in spite that more complex models can be devised 
(for example intermediate phases such as LOFF may be 
present [24l - |26j ). and they would constitute a benchmark 
in testing exotic compact star composition. 

In the rest of this work we quantitatively show that ac- 
curate astrophysical measurements can already constrain 
important parameters of QCD, like the Cooper pair en- 
ergy gap (A) and the strange quark mass {mg), and show 
how this is done given a subset of the present data. 



II. ANALYSIS 

When considering unpaired strange quark matter [uds 
matter) the equation of state (EoS) assumes the simplest 
form when rris — > 0: e = 3P -|- AB, being e, P and B, the 
energy density, pressure and bag constant, respectively. 
The non-zero strange quark mass has the effect of increas- 
ing the energy of the system Q . On the other hand, if 
the system presents pairing interactions between quarks 
it has the net effect of lowering the energy of the system, 
making it more stable. 

There have been different approaches to analyze this 
effect. Here, in order to set a framework for the class of 
models studied, we adopt the MIT bag model thermody- 
namic potential of CFL strange quark matter [1^, [l3 '^^ 



bulk to order A^ , being A the pairing energy gap for the 
phase, given by 
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In this expression, the CFL state is taken as a model 
composed of quark matter without pairing plus a term 
corresponding to the energy of the diquark condensate. 
Being so, the pairing energy gap is taken as a free param- 
eter, instead of being defined in a self-consistent manner 
(see, for example, ,27] for further discussions and more 
consistent models) . The thermodynamic potential for the 
unpaired state is then 



n 



free 



E 



1 

4^ 



n oX O 






i—u,d,s 



(2) 



with fi = (/i„ + /id + fJ's)/S and the common Fermi mo- 
mentum, i/ = 2/X — a/zZ^ _|_ 777,2^2. 

The EoS is finally obtained by relating the energy den- 
sity e to the pressure P as follows 
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where the starred quantities are the dimensionless mass, 
pressure and mass density; x is the dimensionless radius 
and a = ^rrg is an effective parameter which depends on 
the quantities A and ms in a complicated but very mild 
and continuous way. 

Our main goal here was to devise an easy way to dis- 
play the location (locus) of the maximum mass and its 
correspondent radius, e. g., fitting a curve / = f{x,a). 
This locus strongly depends, in our parameterization, 
on a. The higher its value, the higher the maximum 
mass. Because of this, we needed to construct several 
mass-radius relationships varying a within an acceptable 
range, in this case, from 0.20 to 0.36. To comply with 
observations, we needed to reach at least 2Mq. We then 
fitted an expression for the maximum mass {mi,^raax) of 
each sequence and their correspondent radius {xcorr) de- 
pendent on a, yielding the result 



rrii,, 



-O.ieSa^ -I- 0.201a -I- 0.00286 



(8) 



Xcorr = -0.262q;2 + 0.329a + 0.111. (9) 

Restoring the units, we finally find the relation 



with the particle density given by ub — [v^ + 2A'^fi)/n^, 
and the pressure P = —^cfl- 

Due to the fact that iris 7^ 0, the EoS have to be com- 
puted numerically. However, the CFL EoS shows a re- 
markably linear behavior [19|. Because of this, much 
work has been done in order to construct a simple EoS 
that includes the effects of pairing and a heavy strange 
quark, keeping the linear behavior to a certain extent. 
We start our analysis, then, by postulating a linear EoS 
for CFL matter. 

We have parameterized the equation of state for CFL 
matter as 



P = a{c^p - AB.ff) 
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where a is a parameter related with the stiffness of the 
EoS and Be// is certainly related to the bag constant 
B and gap parameter A. Within this approximation, 
it is easy to show that the Tolman-Oppenheimer-Volkoff 
(TOV) and the EoS equations can be written in dimen- 
sionless forms: 
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Here is the point where the dependence on Be// ap- 
pears. Notice that the radius and the mass individually 
depend on Be// but the relation between maximum mass 
and its correspondent radius (when the physical units are 
restored) does not: it depends only on the value of a. 
This feature is already present in Witten's original paper 
and seems to be generic for any linear equation of state. 

From the effective parameters found above, we have 
analyzed the agreement of the mass-radius relation ob- 
tained for different sets of parameters with data for the 
following stars: 4U 1608-522 f20|, EXO 1745-248 [IH, 4U 
1820-30 [13 and PSR J1614-2230 [H]. The results are 
presented in the following section. 
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III. RESULTS 

The value of a as a function of Bg// can be constrained 
by using the observed mass of PSR J1614-2230, consid- 
ering that the maximum mass provided by a specific EoS 



should at least equal the mass of this heavy star. This can 
be obtained from eq. (|5))- ([T(I1) or by varying the values 
of a and -Be// on Eq. (J4|) and constructing the corre- 
sponding mass-radius sequences. Within this previously 
determined region and having generated the mass-radius 
relations for the effective equation of state using these 
constrained values, we compared them with the other 
three stars within a Scr demanded precision. The effec- 
tive parameters, i.e. all the pairs (a, Be//), that would 
be suitable for explaining all four data at the same time, 
i.e. the observed masses and radii and their respective 
error bars, are those corresponding to the line shown on 
the left panel of Fig. [TJ The initial and final points of 
this line were used to generate the mass-radius relations 
shown on the right panel of the same Fig. 

We finally find a very narrow window for a and Be// 
that can provide compatibility with all four measure- 
ments together. As a consequence, the mass-radius re- 
lations shown on the right panel of Fig. [T] are almost 
completely overlapping. 

We would like now to relate the effective parameters 
determined through this analysis with the free parame- 
ters determining the CFL matter state. Unfortunately, 
although the EoS presents a linear behavior, the lin- 
earization of Eq.([T]) is not possible analytically, unless 
some major approximations are made. Without the lat- 
ter, one cannot directly relate a and Be// to A, vtis and 
B. Previous attempts tojind such a relation include the 
work of Alford et al. [23], where a form of parameter- 
ized EoS was proposed, in which QCD corrections to the 
pressure, represented by a parameter c ~ 0.3, have been 
added. This parameterization relied on a power series 
expansion in [i for the thermodynamic potential of the 
quark phase, reading 
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with 04 = 1 — c and ai as free parameters. In their ap- 
proach, the strange quark mass and gap parameter affect 
02 alone and Be// is determined by the density at which 
the transition between quark and nuclear matter phases 
occur, as they were interested in addressing hybrid stars. 
Although they have been able to write ai = 02(7/15, A), 
and one can in fact for CFL quark matter in equilibrium 
(within this framework) separate the thermodynamic po- 
tential expression to order mf with only three terms pro- 
portional to /i^, /i^ and /i°, the chemical potential itself 
has a nonlinear and nontrivial dependence on all other 
parameters, as shown in details in [19]. We have chosen 
not to make this approximation here, so that we are not 
able to write a and Be// as function of rris , B and A in 
analytical closed form. 

Since it is not possible to obtain the expression relating 
the parameterized EoS to the CFL one, in order to verify 
the effectiveness of our own linearization, we have taken 
the full equation of state for CFL strange quark matter, 
as written in Eq. (IT|)-([3]), and analyzed the agreement 



of the mass-radius relation obtained for the same set of 
stars as before. We have covered acceptable ranges for 
the free parameters: the bag constant B > 57 MeV/fm'^, 
the superconductor gap < A < 50 MeV and the strange 
quark mass ttt-s > 100 MeV. Recent measurements give 
a very precise value for the s quark mass at the energy 
scale of 2 GeV: 92.4 ± 1.5 MeV f^g"]. At a smaller energy 
scale this value should be even higher (see [30|, Parti- 
cle Properties - quark masses, and references therein) so 
we decided to adopt the value of 100 MeV as a lower 
boundary. 

After solving the TOV equations for several sets of 
(tos, B, A) and then comparing the values of mass and 
radius with the measurements under consideration, the 
best agreement found is the one giving the mass-radius 
relation shown in Fig. [2] The "optimal" values of the 
parameters found are B ~ 69 MeV/fm'^, A ~ 50 MeV, 
and 7TT-S ^120 MeV (well within the window of stability 
for CFL strange quark matter, thus justifying the quark 
star scenario). Notice, however, that for 4U 1820-30 the 
mass-radius curve is in agreement within '^ 3.4cr only. 
It should also be noticed, nevertheless, that one cannot 
state that in this case the 2>a boundary is the one corre- 
sponding to 99.9% of compatibility, since there are very 
few degrees of freedom involved in the analysis. There- 
fore, the agreement can be considered as acceptable. 

We see that, although the equation of state for CFL 
strange quark matter clearly presents a linear behavior 
in the pressure versus energy density dependency, when 
a parameterization of the type P = aic? p — ABeff) is 
made, the coefficients a and Be// cannot be treated as 
independent from each other. Since we are working with 
a "toy model" , in which one considers a system without 
any pairing plus a term corresponding to the condensate 
energy, A surely adds up to the bag constant but has no 
infiuence on the "unpaired state" . However, the role of 
the strange quar k mass is more subtle. In the analysis 
discussed in [2^, the authors show that Be//, besides 
being related to the vacuum pressure (and, in the CFL 
phase, also to the gap parameter), incorporates an addi- 
tional term which is proportional to m^ when the power 
series expansion is made. So even in a simplified scheme, 
we see that Be// is not dependent on the bag constant 
alone but should also incorporate additional terms which 
are important for determining the stiffness of the EoS. 

In this way, even having found values for the parame- 
ters in the linearization that could fit the data within a 
^ 3cr region, when a full equation of state is employed, 
one cannot reach such a good agreement. This kind of 
problem is also evident in the oldest expression employed 
by Benvenuto and Horvath [3l| ^ = 3^ [c^P - (4 + 6)B] , 
the complicated dependence of the parameters on ttIj, 
B and A (and hence their correlation) affects the accu- 
racy of the Mmax — R locus and blurs to some extent an 
accurate determination of the parameter space. 

Another general aspect worth of noticing is the very 
soft raise of the mass-radius curve for paired strange 
quark matter. It could be attributed to the simplic- 
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FIG. 1: On the left panel, effective parameters allowing compatibility within 3a with recent data (taking into account the 
observed masses and radii and their respective error bars). On the right, mass-radius measurements are shown as crosses with 
3(T error bars. Data from PSR J1614-2230 is shown as a shaded region, since there is only information on the mass for this 
object. The two mass-radius relations shown correspond to the lower and upper limits of a and B^ff that can explain the data 
(shown on the right). 



ity of the model, not considering interactions between 
quarks other than pairing and one-gluon exchange (see 
[?[ for a discussion on the scaling of the bag parameter 
with Qfc ^ 0). A steeper raise of the mass-radius rela- 
tion could indeed provide a better agreement with the 
observed data and perhaps an increase in the maximum 
allowed mass. This could be an important issue since 
it would be impossible to reach masses as high as 2Mq, 
for example, with a single quark matter EoS that would 
also need to explain the values for measured masses and 
radii for the stars 4U 1608-522, EXO 1745-248, and 4U 
1820-30, if they happen ultimately to be accurate and 
correct. 

There have been claims, however, that the values of 
masses and radii for 4U 1608-522, EXO 1745-248 and 4U 
1820-30 have been mistakenly derived ^^. The compu- 
tational method used in the previous analysis was criti- 
cized pointing out possible systematic errors which sub- 
stantially reduce the uncertainties in the masses and radii 
values for these systems. Considering in turn the deter- 
minations given by Steiner, Lattimer and Brown [32|, we 
have repeated the previous calculations to verify whether 
the CFL equation of state is compatible with the recon- 
struction of the data presented by them. The results 
of our analysis for this case are very different, though. 
We see that the very narrow parameter space previously 
found now becomes an extremely wide region of possibil- 
ities, as depicted in Fig[3l 

The error bars for these systems, according to Steiner, 
Lattimer and Brown, are much bigger than previously 
considered and all three stars are clustered with close val- 
ues for their masses and radii (M ~ 1.5M0 and i? ~ 11 
km with more visible differences for the mass when the 
photospheric radius for photon emission is considered to 
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FIG. 2: Best agreement for the mass-radius relation with re- 
cent measurements (data represented as in Fig. [T| consider- 
ing the CFL equation of state. Values for the parameters are 
B ~ 69 MeV/fm^ A ~ 50 MeV, and m^ ~ 120 MeV. 



be much larger than the stellar radius, Rph >> R, the 
color correction fc is a key ingredient for this behavior). 
This leads now to a large set of parameters rris , B and A 
that could in principle explain the data, basically having 
the mass of PSR J1614-2230 as the main quantity pro- 
viding the constrain for the equation of state. In other 
words, one can find several sets of (B, A, nis) that satisfy 
the necessary conditions within the values needed in this 
analysis. The general behavior is to have a lower max- 
imum mass with an increase in the strange quark mass 
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FIG. 3: A second evaluation of the reconstructed data as given in [33] ■ On the left side, mass and radius of 4U 1608-522, EXO 
1745-248 and 4U 1820-30 and mass range for PSR J1614-2230 considering R = Rpu- On the right, the same for R » R^u (see 
text for further details). Two examples of mass-radius relations for CFL stars are shown for the sets of parameters (_B, A, ttIs) 
given by (69 MeV/fm^, 50 MeV, 120 MeV) and (58 MeV/fm^, 62 MeV, 110 MeV) to highlight the fact that, in this case, the 
mass-radius relation is constrained solely by the massive pulsar. 



and bag constant when the values of the other two vari- 
ables are held fixed. As the gap parameter increases, the 
the maximum allowed mass also increases. 



IV. CONCLUSIONS 

Precise measurements for masses and radii of compact 
stars are very important for constraining their composi- 
tion. However, besides being precise, it is necessary that 
they also be accurate for a clear comparison of the mass- 
radius relation coming from very different compositions 
with data. 

Given the present controversy of the radius determi- 
nation quoted in the paper, this appears to be the main 
issue to be resolved for an advance of the conclusions of 
studies like ours. Meantime we can only conclude that 
only with small error bars for both mass and radii deter- 
minations (as in the analysis of [20l - |22| ) the parameters of 
the CFL equation of state can be effectively constrained. 
Otherwise, several sets of parameters defining the CFL 
EoS can be used to match the data. In this case, more 
advanced detection techniques and physical analysis are 
required to constrain the radius Rph in X-ray burst mod- 
els. 

Since one should not treat a and -Be// as indepen- 
dent quantities in a parameterized EoS for CFL strange 



quark matter like the one proposed here, an analysis 
performed with a parameterization with correlated co- 
efficients should be explored. We intend to address this 
subject in a future work. 

An additional interesting report is the very low mass 
of at least one neutron star mass observed (M = 0.87 ± 
O.O7M0 for 4U 1538-52) (sl, actually not well explained 
by ordinary stellar evolution arguments. This "low-end" 
mass extreme is also important for the quest of exotic 
compositions in their interiors. In fact, confirming such 
a measurement would imply a radius of i? < 9 km within 
the self-bound hypothesis, while the measurement of, say, 
i? > 11 km would be totally inconsistent with it, pointing 
to a normal composition. We are at the edge of finally ad- 
dressing the fundamental composition of compact stars, 
at least on a first level that leaves a handful of possibili- 
ties left out of a huge number of proposals. 
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